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Abstract 



We show that the decay of a passive scalar 9 advected by a random incom- 
pressible flow with zero correlation time in the Batchelor limit can be mapped 
exactly to a certain quantum-mechanical system with a finite number of de- 
grees of freedom. The Schrodinger equation is derived and its solution is an- 
alyzed for the case where, at the beginning, the scalar has Gaussian statistics 
with correlation function of the form e~l^~^l . Any equal-time correlation 
function of the scalar can be expressed via the solution to the Schrodinger 
equation in a closed algebraic form. We find that the scalar is intermittent 
during its decay and the average of (assuming zero mean value of 0) falls 
as e"'*'"^* at large t, where D \s a. parameter of the flow, 7^ = \ct{'d — a) for 
< a < 3, and 7o = | for a > 3, independent of a. 
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Kolmogorov theory (K41) remains the cornerstone of our understanding of fully de- 
veloped turbulence. This simple theory predicts a scaling law (the famous Kolmogorov- 
Obukhov k~^^^ law) of the energy spectrum that is in remarkable agreement with experi- 
mental data. Since the 1980 's, however, data gathered have consistently pointed out the 
failure of K41 in predicting the scaling law of high-order correlation functions [0,0. The 
breakdown of K41 is closely related to the non-Gaussianity of the distribution of velocity 
increments. The phenomenon, dubbed intermittency, has become one of the central issues 
of theoretical works on turbulence. Recently, it has been found that the intermittency of 
a passive scalar advected by a turbulent flow might be even stronger than that for the ve- 
locity [^. Such observations have led to the hope that the study of simple models, such as 
the Kraichnan model of scalar advection (see Refs. |5HT0| and below), may provide clues to 



understand the much more complex Navier-Stokes intermittency. 

In this paper, we consider the problem of turbulent decay of a passive scalar. In other 
words, we want to find statistical properties of a scalar satisfying the equation 

dtO + VidiO = kAO, (1) 

where k is a small diffusivity, Vi is a Gaussian random field, which is white in time, 

(t;,(t,x)t;,(t',y))=5(t-t')/.,(r), (2) 

and 

^■(r) = V6,, - D (^5.,r« - r^'W.^j , (3) 

where r = x— y, and C, is some real number. The Kraichnan model usually contains a random 
external scalar source in the right-hand side (RHS) of Eq. (|l]). Such a source would make the 
steady state possible, but since we are interested in the decay, it is assumed that the source 
is absent. We will, furthermore, turn our attention to the Batchelor limit = 2, which 
corresponds to smooth flows with very large velocity correlation lengths (for comparison, 
the inertial range of real turbulence corresponds to = |.) This limit has attracted recent 
interest due to its good analytical features ||TT| , p!2[| . 

Our result is that the scalar becomes more and more intermittent during the decay. 
Specifically, we found that the average of (|^(x)|"), where a is an arbitrary positive number, 
decays as e'^"^^ at asymptotically large t, where 7^ = |a(6 — a) if a < 3 and 7^ = | when 
a > 3. The flatness {9^) / {6"^) ~ ^^^1^ goes to 00 as t grows. This is in sharp contrast with 
the steady-state case, where the scalar statistics is largely Gaussian Wl . 



To attack the problem, we will reduce it to a certain problem of quantum mechanics, 
which can then be solved (for another attempt to apply quantum mechanics to turbulence, 
see [|13|.) We first note that the probability distribution functional of the scalar, which will 



be denoted \E'[t, ^], can be expressed in term of a path integral ||l4| 



^\t,Q\ = JvTT{t,x.)Ve{t,x.)Vvi{t,x.) p[v] exp i j dt d^T^^dtO + VidiO - kAO) 



(4) 



where the Gaussian measure for the velocity p[v\ is chosen to satisfy Eq. (0). The auxiliary 
variable tt enforces Eq. ([1|). Integrating over v, one obtains 
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'^{t,e) = jvnVe exp i jdxTidte-]^jdtdxdyT^{t,x)die{t,x)fij((K-y)Ti{t,y)dje{t,y) 
—IK / dx7cA9 . 



The path integral describes the evolution in Euclidean time of a quantum field theory with 
the Hamiltonian [|l^] 

H = ^J rfxt/y 7r(x)a,^(x)^ (x - y)vr(y)9,5(y) +tf^J dxirAO, (5) 

where 6 and vr are conjugate variables satisfying the usual commutation relation [^(x), 7r(y)] = 
i5(x — y). The operator ordering in Eq. (^) corresponds to the physical regularization of 
the path integral @. The evolution of the distribution functional "^[6] is described by the 
Euclidean version of the Schrodinger equation, dt"^ = —H"^. Note that the functional \E' 
itself, not its square, determines the probability distribution of 6. The average of, e.g., |6'|" is 
defined as (|^|") = fDO |^|"\E'[^]. In further discussion, we will use the quantum-mechanical 
terminology, so the terms "probability distribution functional" (PDF) and "wave function" 
are used interchangeably. 

In the Batchelor limit (|^), the Hamiltonian can be simplified considerably. We will 
concentrate our attention to the homogeneous case, i.e., when the system is invariant under 
spatial translations. In the quantum language, this means that we restrict ourselves to the 



states |\I') having zero total momentum, Pi\'^) = 0, where Pi = /(ix7r(x)(9j6'(x) [0. With 
this restriction, the Hamiltonian (^) can be rewritten into the following form: 

H = —(^ALijLij LiiLjj LijLji) + ikDh (6) 

where the operators Lij and Dij are defined as 

Lij = J dxxiTT{x.)dj9{'x.), Dij = j d'x.T:{i()didj6iyi). 

It is straightforward to check that Lij and Dij form a closed algebra with the commutation 
relations, 

[Lij, Lki] = i{SjkLii — SuLkj), 

[Lij, Dki] = -i{5iiDjk + SikDji), (7) 
[Dij,Dki] =0. 

The fact that the algebra is closed implies that the system is actually one with a finite number 
of degrees of freedom. The quantum field theory thus degenerates to quantum mechanics. 
Notice that Lij form a closed subalgebra. Indeed, they are the operators of linear coordinate 
transformations. In fact, only the SL(3,R) generators enter the Hamiltonian (^) (cf. [0.) 
H is invariant under the SO (3) algebra formed by the antisymmetric part of Lij. 

In principle, the Schrodinger equation with H defined in Eq. can be solved (at least 
numerically.) In this paper, we will choose a representation of the algebra (|^) where H has a 
relatively simple form, but the physics is nontrivial. Our choice is inspired by the observation 
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by Townsend fl^ that a Gaussian shape hot spot preserves its Gaussianity when advected by 



Batchelor-hmit velocity flow (for a somewhat similar discussion without quantum mechanics, 
see |]18[.) Let us for a moment concentrate on the states in which 9 has Gaussian statistics. 
This corresponds to the wave functions of the form '^[9] ~ e'x.Y>{—\9K~^9), where K{x — y) = 
{9{x)9{y)). We will further restrict ourselves on functions K that have the Gaussian shape, 
K{x — y) ~ exp[—jbij{x — y)i{x — y)j]. More strictly, we require that, in Fourier components, 
the spectrum of 9 has the form 



'ik)9ik')) = 9oexp(^~a,,hk,>jS{k - k'). 



where is a constant independent of aij = ipij) ^ (one can choose 6*0 = 1.) Denote such 
states as \aij). The group elements act on \aij) as follows: 

e-^-^-|a,,) = |e-'^a(e-T) if A. = 

= |ai,- + 4A,-). (8) 

We now choose our representation to be the one acting on the subspace of the Hilbert 
space that contains all linear combinations of \aij) (although the latter do not form an 
orthogonal basis.) A vector in this subspace is characterized by the function iplaij), which 
is the coefficient of the expansion |\&) = Jdaij'ilj{aij)\aij). In general, the scalar statistics 
in is not Gaussian. The operators Lij and Dij can be written as first-order differential 
operators with respect to aij, and the Schrodinger equation becomes a second-order PDE 
on ip. 

Moreover, if the initial condition is isotropic, i.e., invariant under S0(3) rotations eijkLj^, 
the wave function depends only on the eigenvalues of the matrix aij, not on the Eulerian 
angles characterizing the orientation of the eigenvectors. The wave function is now a function 
of three variables, 'ip{ui,U2,U3), where we have denoted the eigenvalues of aij as e^"\ We 
re-scale ip so that the state is expressed via ip{u) as 

\^!) = jduidUi){u)\a{u,U)), (9) 

where a{u,U) = t/diag(e^"')[/~^, U belongs to S0(3), and the integration over U is per- 
formed using the invariant measure on the SO (3) group manifold. 

The Schrodinger equation ip{u) can then be derived (details are found in [0). It has 
the form 

dt^ = D{dl + dl + dl - d,d, - d2ds - d^d,)^!; - [3Dd,{m + 2K5,(e-2«'^)], (10) 

i=l 



where di = d/dui 



g4ltl _ g2(M2+M3) 

fl = f{Ui;U2,U3) 



f2 = fiu2;us,ui), h = fiu3;ui,U2). (11) 

Special caution is required when two of Ui are equal to each other, however this will not 
affect our subsequent discussion. 
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To fully define the problem, the initial condition of ip{u) is needed. One can take as the 
initial state the vector \aij), where Ojj = diag(l, 1, 1). This corresponds to a scalar that has 
Gaussian statistics, zero mean value, and the correlation function {6{x)6{0)) proportional 
to at t = 0. The correlation length of 9 is taken to be of order 1. In terms of the 

initial condition is ip{t = 0,u) = 6{ui)6{u2)S{u3). 



Equation (|l^) can be interpreted in an intuitive way by using a three-dimensional random 



walk that has the Fokker-Planck equation coinciding with Eq. (|TD|) | |2Cm , 

Mi = 3Dfi + 2Ke-^^^ + (12) 
where are white noises that correlate as follows: 

6 + 6 + ^3 = 0, 

{mut')) = iumt')) = imut')) = '^mt - f). (13) 
{mut')) = m)ut')) = m)ut')) = -D6{t-t'). 



Let us discuss the physical meaning of Eq. (|T2|). A point {ui,U2,u-i) corresponds to the 
configuration of 9 having the spectrum (|6'(k)p) ~ exp(— | ). In the configuration 

space, 9 is approximately constant inside an ellipsoid with major axes proportional to 
When advected by the flow, this ellipsoid is subjected to random linear transformations. If 
the only transformations of the ellipsoids are those that stretch or compress the ellipsoid in 
the directions of its major axes, the results would be iii = ^i, where are random. Equa- 
tion ([13D reflects the conservation of the volume of the ellipsoid during random stretching 
and compressing. However, the ellipsoid may be subjected to stretching or compressing in 
directions other than the major axes, as well as to shearing. These effects are accounted for 
by the term 3Dfi on the RHS of Eq. ([T^). The incompressibility is not violated, due to the 
identity /i + /2 + /s = 0. The terms 2Ke~^"' are not important unless one major axis of the 
ellipsoid is as small as the diffusion scale. In the latter case, diffusion smears out the scalar 
and causes it to be correlated at a larger distance. This is exactly the effect of the 2fi:e~^"' 
terms in the Langevin equation. Due to the sign of these terms, the volume of the ellipsoid 
and, hence, also Ui + U2 + u^, always grows during the random walk. 

Since any correlation function can be computed for \aij), where the scalar statistics is 
Gaussian, one can find any correlation function with respect to if one knows the solution 
to Eq. (|T0|) (e.g., from numerical integration.) For example, the average of |^|" (a > 0) over 
the state \a{u,U)) is proportional to e~"("i+"2+«3)/2. therefore, its average with respect to 

I*) is, 



(l^r) = C^{9\t = 0)r' / dui^iu) exp 



--(Ml + M2 + U-i) 

where C„ = 'K-^l'^2'^l'^V{{a + l)/2]. This relation is exact. 

When K is small, the exponential behavior of (|6'|") can be found analytically. This can 
be done by using the path-integral description of the random walk (O) and finding the 



saddle-point trajectories that dominate \9\^ [|l^]. In this paper, we use a heuristic, yet more 



physical, method to find the large time behavior of 

Let us assume that after letting the system (|12D evolve for a while, the values of mi, M2, 
and M3 become widely separated. We assume U\ < < M3, and wide separation means 
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U2 — Ui ^ 1, M3 — M2 ^ 1. From Eq. (0) one sees immediately that in this regime /i = —1, 
/2 = 0, and /s = 1 (in fact, these asymptotic values of are related to the Lyapunov 
exponents, see, e.g., Ref. pT| ). 

Let us first ignore the term proportional to diffusivity in Eq. (p!2D. The velocity iii has 
two contributions: one from /j and another from the noise ^j. The first contribution implies 
that the mean values of Ui drift with constant velocities, ui{t) = —3Dt, U2{t) = 0, and 
^3(^) = 3Dt, while the noises make Ui fluctuate around these mean values. The condition of 
wide separation of m's is satisfied when t ^ D~^. The advection, on average, compresses a 
fluid element in one direction by a factor of e^^* and stretches it in another direction by the 
same factor. The remaining third direction is not substantially compressed or stretched. In 
this regime, the diffusion is still not operative, and (|^|") remains constant. 

At t = {6D)~^\n D^^ if k is very small), the mean value of ui becomes ^Ihk. 

The term /te~^"^ in the Langevin equation (^) cannot be ignored anymore. Physically, 
regions of different 6 have been brought this close together so that diffusion is no longer 
negligible. Let us consider the equation for ui, ui = —3D + 2/te~^"^ +^1, 

k In K. 

The first term on the RHS pushes ui toward smaller values, while the second term prevents 
Ml from becoming substantially smaller than Wmin- The variable ui thus fluctuates around 
Mjnin- Therefore, the random walk becomes effectively two-dimensional: 

imUt')) = imUt')) = 2D5it - t'), (14) 
{Ut)Ut')) = -D6it-t'). 

Additionally, it is required that U2 + not decrease with time, due to the previously 
found fact that Ui + U2 + M3 can only increase (if U2 + Ms decreases, this means that Ui 
steps away from the value Ui = Umin-) Now there is a possibility for 16*1" to decay, since it 
is proportional to e~"("i+"2+«3)/2^ ]^^^ _l_ ^2 + is no longer a constant. Assuming that 
the random walk ([T^) starts at U2 = U2 and = M3, the distribution of U2 and M3 at large 
times is Gaussian: 

p{u2,us) ~ exp 



3Dt 



{U2 - Ma)^ + ("3 - - ^Dtf + {U2 - M2)(W3 " ^3 " ^Dt) 



(15) 



The mean value of 1^1°" can be computed by taking the average of e~"("2+M3)/2 ^-y^j. ^j^g distri- 
bution (pISl). Consider the case of < a < 3 first. The integral Jdu2 du^ p{u2, M3)e~"'-"2"'""^)/^ 
is dominated by the region near M2 — M2 = —\aDt, u^ — u^ = (3 — a/2)Dt. The value of the 
average is proportional to e"'^"-^*, where 7q, = \a{6 — a). 

Note that the region where the integral is saturated has U2 decreasing with time, U2 = 
^aDt. Eventually, U2 will become as small as Mmm; and the term Ke"^"^ in the Langevin 



u 



equation becomes important. Now, both ui and U2 fluctuate around Mmin- However, as we 
will explain, the exponential decay law does not change. Indeed, when ui and U2 remain 
approximately constant, the evolution of U3 is described by the one-dimensional random 
walk, 

us = SD + ^3, {Ut)Ut')) = '^D5{t - t')- 
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The distribution of U3 is now p{u3) ~ exp[— (4Dt)~^(-U3 — -Ug — 3Dt)'^]. Taking the average 
of e~""^/^ (which is proportional to (|^^|") since Ui and U2 are constant), one finds that the 
decay law is still e~'^"^^ where 7^ = \a{Q — a). 

For the particular case a = 2, our result can be checked against the calculations based 
on the exact evolution equation for the scalar spectrum This comparison has been done; 
the results indeed agree. 

When a > 3, the solution U2 ~ —^aDt, tis ~ (3 — a/2)Dt is no longer realizable, 
since it has decreasing U2 + u^. The average of |^|" is then determined by the edge of the 
distribution function, i.e., by U2 ~ and ~ f-Di, or, after U2 reaches Umin, U2 ~ Mmin 

and U3 ~ const. The expectation value decays as e"^-^*/^. The reason the decay law does 
not contain a is the following: when a > 3, the main contribution to comes from the 
realizations in the statistical ensemble where 6 is unaffected by diffusion (i.e., the ellipsoid 
in which 6 is approximately constant has never been too thin during its evolution.) The 
average (|^|") is thus determined by the probability of such realizations, which depends only 
on characteristics of the flow but not on a. This probability, as has been found, falls as 
g-9Dt/4^ This implies, in particular, that the flatness (9^) / {9"^) grows as e'^^*/^, meaning that 
the scalar becomes more and more intermittent during its decay. 

More careful analysis shows that the decay law e~"'°'^^ that we have found is valid only 
at large enough t. At intermediate t, there is a smooth transition from (|^|") = const to 
d^^l") ~ Q-iaDt |j|g|^ T^liQ full analysis does not change the long-time tail of (|^|"). 



In conclusion, we have shown that by mapping to quantum mechanics, the problem of 
turbulent decay of a randomly advected scalar in the Batchelor limit can be made completely 
solvable. The power of the approach described in this paper is not limited to the calculations 
of (|6'|°); analogous calculations can be done for any equal-time correlation function. For 
example, the long-time tail of {\dx9\°') is also e~^"^^ with the same 7q,. The situation here 
is not similar to the steady state, where the scalar and its derivatives have very different 
statistics, with the scalar being largely Gaussian and its derivatives being intermittent W2 



The relevance of the techniques presented and results to the general problem of intermittency 
is yet to be explored. 

The author thanks E. Farhi, J. Goldstone, and K. Rajagopal for helpful discussions, and 
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